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ABSTRACT 


An  analytic  and  experimental  program  was  undertaken  to  define 
the  near  wake  structure  behind  flapped  and  unflapped  wings.  The 
vortex  wake  structure  is  determined,  given  the  wing  lift  and  drag 
distributions  with  models  developed  in  the  spirit  of  Betz.  A 
procedure  to  estimate  the  turbulent  kinetic  energy  distribution  in 
the  vortex  is  also  given,  as  well  as  a method  to  determine  discrete 
vortex  positions  in  the  downstream  wake.  Theoretical  models  are 
compared  with  detailed  three-component  velocity  measurements  in  the 
wake  of  a flapped  model  wing.  Predicted  wake  velocities  are  in 
very  good  agreement  with  measurements  made  in  the  wake  at  distances 
d-'rfnstream  where  roll-up  is  essentially  complete.  Computed  dis- 
crete vortex  positions  in  the  downstream  wake  compare  most  favor- 
ably with  results  of  flow  visualization  studies.  Estimates  of  the 
downstream  turbulent  structure  of  the  wake  vortices  are  made  using 
an  invariant  turbulent  model  developed  at  A.R.A.P. 
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SECTION  I 
INTRODUCTION 


The  hazard  associated  with  aircraft  wake  turbulence  is 
now  well  known,  and  there  is  currently  underway  an  extensive 
effort  by  the  PAA  and  NASA  to  respond  to  this  problem.  This 
effort,  however,  is  not  concerned  directly  with  vortex  wake 
problems  which  are  unique  to  the  Air  Force;  namely,  hazards 
associated  with  short  Interval  take-offs,  mid-air  refueling, 
and  formation  flight.  Since  wake  intensity  can  be  so  severe 
that  encountering  aircraft  have  been  known  to  lose  control,  it 
is  of  grave  importance  to  be  able  to  determine  for  which  air- 
craft and  under  what  operating  conditions  hazards  to  other 
aircraft  exist.  In  regard  to  the  specific  operating  conditions 
which  are  unique  to  the  Air  Force,  this  hazard  can  only  be 
assessed  from  a detailed  description  of  the  aircraft  wake.  This 
description  must  answer  three  basic  questions: 

1)  What  is  the  pattern  of  vorticity  shed  in  the  Immediate 
vicinity  of  an  aircraft,  and  how  does  this  vorticity 
pattern  tend  to  concentrate  (roll  up)  behind  the 
aircraft? 

2)  What  are  the  effects  of  diffusion,  turbulent  and 
laminar,  on  the  vortex  patterns  that  develop?  We 
call  these  effects  aging. 

3)  Are  these  initial  patterns  of  concentration  stable, 
or  will  the  patterns  undergo  Crow  instability  and 
vortex  breakdown? 

The  answers  to  all  of  these  questions  are  relevant  to  the  Air 
Force . 


The  research  described  in  this  report  attempts  to  provide 
a logical  answer  to  the  first  two  questions.  In  Section  II,  we 
review  the  Betz  roll-up  model  and  Its  extension  to  include  the 
roll-up  of  '’interior"  vortices.  While  these  techniques  are  now 
well  documented  in  the  literature,  we  include  them  here  and  show 
how  a new  physical  interpretation  of  the  model  allows  further 
generalization.  The  result  is  a model  which  can  predict  the 
vortex  axial  and  swirl  velocity  distribution  for  either  tip  or 
"interior"  roll-ups  given  the  wing  lift  and  drag  distribution. 

The  structure  of  each  individual  vortex  in  the  wake  Is  not  com- 
plete, however,  without  specifying  how  the  turbulent  kinetic 
energy  in  the  boundary  layer  on  the  wing  might  be  distributed  in 
the  rolled-up  vortices.  A procedure  for  doing  this,  in  the  spirit 
of  the  Betz  model,  is  also  presented.  The  local  turbulent  struc- 
ture of  the  individual  vortices  may  then  be  calculated  using  a 
three-dimensional  turbulent  vortex  program  which  has  been  devel- 
oped by  A.R.A.P.  for  ARL  under  Contract  F33615-72-C-2116 . The 
structure  of  the  wake  is  completed  by  developing  a model  by  which 
the  motion  of  the  vortices  which  comprise  the  wake  may  be 
determined.  In  short.  Section  II  develops  theoretical  models 


from  which  a complete  description  of  the  near-field  rolled-up  air- 
craft vortex  wake  is  possible. 

Section  III  describes  the  design  and  implementation  of  a 
test  program  to  check,  where  possible,  the  above-described  models. 
The  test  program  was  carried  out  in  the  V/STOL  tunnel  at  NASA’s 
Langley  Research  Center.  Detailed  measurements  were  made  by 
Langley  personnel  of  the  wake  structure  and  the  model  wing  lift 
distributions . 

In  Section  IV,  the  results  from  this  test  program  are 
analyzed  and  compared  with  calculations  from  theoretical  models. 

Finally,  in  Section  V,  conclusions  and  recommendations  are 
offered . 
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SECTION  II 

REVIEW  AND  DEVELOPMENT  OF  THEORETICAL  MODELS 


It  is  well  known  that  a finite  aspect  ratio  wing  which 
develops  circulation  lift,  sheds  vorticity  as  a consequence  of 
this  lift.  This  vorticity,  which  forms  the  wake,  is  shed  from  the 
trailing  edge  and  is  commonly  referred  to  as  a shed  vortex,  sheet 
(the  fluid  containing  vorticity  is  very  thin  in  a direction  normal 
and  aft  relative  to  the  planform  of  the  wing).  At  the  trailing 
edge,  the  sheet  is  more  or  less  planar,  but  does  not  remain  so. 

A short  distance  downstream,  the  sheet,  as  a result  of  a convec- 
tive instability,  tends  to  roll  up  into  discrete  vortices,  as 
shown  schematically  in  Figure  1.  These  vortices  can  be  quite 
persistent  and,  when  of  sufficient  intensity,  are  a hazard  to 
other  aircraft . 

Calculation  of  the  details  of  the  roll-up  is  quite  difficult, 
and  investigators  have  resorted  to  simplified  models  to  describe 
the  phenomenon.  One  such  model  popularized  by  Westwater  in  1935 
(Rof.  1)  calculates  the  roll-up  of  a two-dimensional  sheet  by 
replacing  the  sheet  with  discrete  vortex  elements.  This  approach 
has  formed  the  basis  of  several  recent  investigations  (Refs.  2-5). 
One  difficulty  is  that  the  similarity  solution  of  Kaden  (Ref.  6) 
shows  that  the  center  of  the  rolled-up  spiral  contains  an  infinite 
number  of  turns  and,  as  Westwater  points  out,  can  never  be 
modeled  by  a finite  number  of  vortices.  Recently,  Moore  (Ref.  5) 
has  used  a scheme  whereby  the  difficulties  associated  with  model- 
ing the  spiral  structure  with  discrete  vortices  are  circumvented. 

He  replaces  the  exact  spiral  structure  with  an  irrotational  tip 
vortex,  thereby  eliminating  the  need  to  keep  track  of  the  many 
discrete  vortices  which  model  the  spiral.  While  the  results  of 
these  calculations  are  encouraging,  it  is  unfortunate  that  the 
detailed  structure  Is  lost.  Calculations  based  on  this  method  are 
valuable  in  that  they  are  able  to  estimate  time  to  roll-up. 

When  the  details  of  the  roll-up  are  not  needed,  two  models 
have  been  suggested  to  obtain  the  vortex  wake  structure.  The 
first  model,  proposed  by  Spreiter  and  Sacks  (Ref.  7),  equates  the 
3wirl  kinetic  energy  per  unit  length  of  wake  to  the  Induced  drag 
of  the  aircraft.  The  calculation  requires  an  assumption  as  to  the 
nature  of  the  swirling  velocity  distribution,  with  sufficient  free 
parameters  that  circulation  about  each  vortex  and  the  impulse  of 
the  system  are  conserved.  In  Reference  7,  the  vorticity  was 
assumed  to  be  distributed  uniformly  in  the  vortex.  For  an  ellip- 
tieally  loaded  wing,  vortex  radius  was  obtained  to  be  0.155  the 
semi-span  of  the  wing.  While  calculations  of  this  nature  are 
straightforward,  they  do  not  give  a unique  relationship  between 
wing  load  distribution  and  vortex  velocity  distribution. 

The  second  model  was  proposed  by  Betz  and,  while  available 
for  some  time,  received  little  attention  until  Donaldson  (Ref.  8) 
showed  that  the  swirl  velocity  distribution  calculated  In  this  way 
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Figure  l.  Illustration  of  the  roll-up  of  shed  vorticity  from  a simply  loaded  wing 
(from  Ref.  8) 
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compared  most  favorably  with  measurements.  Recently,  several 
studies  in  the  spirit  of  Betz  have  been  undertaken  (Refs.  9-12), 
and  Donaldson,  et  al.  (Ref.  13)  has  shown  how  the  roll-up  of  flap 
vortices  may  be  calculated  according  to  the  Betz  assumption. 

Since  the  comparison  of  the  Betz  model  and  the  extensions  proposed 
in  Ref.  13  with  experimental  measurements  forms  a significant 
portion  of  the  effort  reported  herein,  it  is  appropriate  to  review 
this  model. 
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1.  THE  BETZ  ROLL-UP  MODEL  AND  EXTENSIONS  TO  INCLUDE  INTERIOR 
VORTICES 


The  method  described  by  Betz  for  calculating  roll-up  relates 
the  circulation  r at  wing  station  y to  the  circulation  I” 
calculated  at  radius  r in  an  axisymmetric  line  vortex.  The 
method  is  based  upon  the  assumption  that  global  invariants  appli- 
cable to  an  unbounded,  two-dimensional,  incompressible,  inviscid 
fluid  medium  may  be  applied  locally  behind  a wing  to  obtain  an 
approximate  description  of  the  vortical  wake.  The  fundamental 
assumption  Is  that  vortical  motions  are  such  as  to  preserve  the 
inertial  moment  about  the  centroid  of  vorticity;  that  is. 
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where  y(y) 
y and  s . 

is  the  centroid  of 

the  vorticity  shed  between 

stations 

Equation  (1)  is  approximate  and,  as  will  be  shown  in  Section 
II. 2,  can  be  manipulated  so  as  to  allow  physical  Interpretation. 
With  Eq.  (1)  and  a statement  of  Kelvin’s  theorem 


f’(r) 


(3) 


which  is  exact  for  an  inviscid  flow,  Betz  was  able  to  give  a 
rather  complicated  expression  for  the  swirling  velocity  in  the 
rolled-up  wake  of  an  elliptically  loaded  wing.  By  manipulating 
Eqs.  (1)  - (3),  Donaldson,  et  al.,  Rossow,  and  Jordan  have  inde- 
pendently shown  the  surprisingly  simple  result:  the  relationship 
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between  r an  <2  y is 


r a y(y)  - y (*0 

This  result,  along  with  Eq.  (3),  states  that  the  value  of  the 
circulation  at  wing  station  y is  the  value  of  the  circulation 
at  radial  distance  r in  an  axisyrametric  vortex.  The  radial 
distance  r is  equal  to  the  distance  from  y to  the  centroid  y 
of  all  the  shed  vorticity  outboard  of  y . When  all  the  vorticity 
can  be  considered  rolled  up,  the  vortex  center  is  located  at 
y = y(0)  In  order  to  preserve  the  vex*tical  momentum  of  the  flow. 
Since  r also  equals  y(Q)  , the  circular  regions  containing 
vorticity  Just  touch  along  the  aircraft  centerline.  Figure  2 
depicts  the  Betz  roll-up  model. 

When  an  aircraft  has  flaps  and/or  spoilers  deployed,  the 
wing  load  distribution  may  be  such  that  a single  roll-up  proceed- 
ing from  the  tip  is  no  longer  possible.  In  such  cases.  Ref.  13 
has  given  a criterion  which  specifies  how  the  vorticity  distribu- 
tion will  divide  itself  and  roll  up  into  multiple  discrete 
vortices.  The  criterion  is  best  presented  by  considering  the 
example  given  in  Figure  3.  The  function  dr/dy  is  the  strength 
of  the  vortex  sheet  shed  from  the  wing.  Consider  the  distribution 
of  the  absolute  value  of  the  shed  vorticity  j dr/dy \ . This 
function  has  three  minima  which  are  denoted  by  points  A,  B,  and  C. 
It  was  assumed  that  the  vorticity  shed  between  B and  C would 
roll  up  Into  what  was  generally  called  a fuselage  vortex.  The 
vorticity  shed  between  stations  A and  B would  form  a flap 
vortex,  and  all  the  vorticity  outboard  of  A would  roll  up  ar  a 
tip  vortex.  The  strengths  of  these  vortices  are  rc  - rg  , 

fB  “ * and  » respectively.  Furthermore,  it  was  assumed 

that  the  roll-up  of  the  flap  and  fuselage  or  "interior"  vortices 
would  proceed  from  the  points  where  1 dr/dy | is  maximum,  namely, 
ytnR  and  ymC  . This  Is  physically  motivated  and  compatible  with 

observation.  Recently,  Yates  (Ref.  14)  has  checked  these  assump- 
tions by  the  direct  calculation  of  the  initial  in-plane  accelera- 
tion of  the  vortex  sheet.  He  was  able  to  show  that  these 
assumptions  are  quite  accurate  and  give  a simple  method  to  deter- 
mine the  number  and  strength  of  the  vortices  to  be  expected  in  a 
>rake. 


The  circulation  distribution  in  "interior"  vortices  was 
calculated  In  Ref.  13  by  extending  the  Betz  model.  To  see  how 
this  was  done,  consider  the  vorticity  shed  between  stations  A 
and  B in  Figures  3 and  4.  First,  the  assumption  regarding  the 
relationship  between  the  inertia  moment  of  the  vorticity  distri- 
bution (Eq.  (1))  i3  modified  to  read 
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Figure  2,  Illustration  of  the  Betz  roll-up  model  for  a simply 
loaded  wing  (from  Ref.  33) 
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where  and  y2  are  arbitrary  points  inboard  and  outboard  of 

ymB  * y12  ls  deflned  by 


(6) 


and  is  the  centroid  of  the  shed  vortlcity  between  y^  and  y2  . 
Kelvin’s  theorem  rewritten  for  an  interior  vortex  is 


r(yx)  - r(y2) 


r 

Jo 


d; 


r(r) 


(7) 


Equations  (5)  - (7)  are  not  sufficient  to  describe  an  ’’interior” 
vortex  until  a relationship  between  y^  and  y2  is  prescribed. 
This  relationship  was  taken  to  be 

(y2  “ y12)2  = *yl  “ y12)2  <8) 

and  is  equivalent  to  the  expression  found  by  Betz  for  a tip  roll- 
up. By  manipulating  Eqs.  (5)  - (8),  it  may  now  be  shown  that 


r - y2  - yi2  - y12  - yj  (9 

The  similarity  with  the  result  obtained  for  tip  roll-up  (^)  is 
obvious.  Equation  (9)  taken  with  (7)  determines  the  circulation 
distribution  in  an  "interior"  vorte?;,  providing  y^  remains 
greater  than  y™  and  y2  remains  less  than  yA  . If  y2 
reaches  yA  before  y^  reaches  , the  roll-up  is  to  be  con- 
tinued by  holding  y2  and^  hence,  r2  fixed  and  proceeding  with 
the  roll-up,  letting  r = y^2  - y^  until  y^  reaches  yg  . If 
yi  reaches  yg  before  y2  reaches  yA  , then  roll-up  is 
continued,  holding  y^  and,  hence,  fixed,  letting  r * y2  - 

y12  * 


Before  proceeding  to  the  next  section  and  showing  how  axial 
velocity  in  the  vortex  might  be  included  in  the  Betz  model,  we 
will  present  the  results  obtained  in  Ref.  13  regarding  the  magni- 
tude of  the  swirling  velocity  at  the  center  of  the  vortex.  It 
was  shown  that 


i 
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V(0)  * - - (10) 

where  dr/dy  is  to  be  evaluated  at  y * b/2  for  a tip  vortex 
and  at  (the  point  of  maximum  sheet  strength)  for  an  "interior" 
vortex.  This  result  ls  significant  in  that  the  maximum  inviacid 
swirling  velocity  to  be  expected  is  simply  a function  of  the 
maximum  rate  of  change  of  the  load  distribution.  Of  course,  these 
velocities  are  never  achieved  in  a real  vortex  since  viscous 
effects  require  that  V(0)  * 0 . 
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2.  THE  EFFECTS  OF  WING  DRAG  ON  INVJSCID  VORTEX  STRUCTURE 

While  the  Bet*  model  and  extensions  to  Include  ’’Interior* 
vortices  provide  much  information  regarding  the  mean  discrete 
vortex  structure , it  is  not  complete  until  we  determine  how  the 
wing  drag  distribution  influences  this  structure.  The  further 
extension  of  the  Betz  model  to  include  this  effect  is  straight- 
forward  and  can  be  physically  motivated.  The  development  can  be 
shortened  by  considering  only  the  roll-up  of  an  "interior"  vortex. 
The  result  for  a tip  vortex  is  obtained  by  setting  the  upper  limit 
of  integration  y2  equal  to  b/2  in  the  integrals  which  follow. 

Equation  (5)  is  multiplied  by  -pU^/2  and  integrated  by 
parts  to  yield 


pUm  f 2 

~2~  r2<y2  " y12J 


“ ri^yl  " y12*  + rBr^J 
ry  2 

♦ / i(o)(n  - y12)dn 
Jyi 

* J*  PCVtOU^irc  d ; 
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where  l(y)  « -pU  I*(y)  is  the  sectional  wing  loading  exerted  on 
the  fluid.  The  first  term  in  Eq.  (11)  vanishes  when  (7)  and  G) 
are  substituted.  The  remaining  terms  prescribe  the  distribution 
of  angular  momentum  in  the  vortex.  The  Betz  model  therefore 
distributes  the  angular  momentum  such  that  the  torque  exerted  by 
the  wing  (calculated  about  y12)  between  yi  and  y2  equals  the 
flux  of  angular  momentum  through  a circle  of  radius  r . In  light 
of  this  physical  Interpretation,  nonunlform  axial  velocity  in  the 
vortex  may  be  included  by  modifying  Eq.  (11)  to  read 


i(n)(n 


y12)db 


PtVU)UU)2*c  d? 


(12) 


The  geometry  of  the  flow  model  for  a tip  roll-up  might  be  thought 
to  be  that  illustrated  in  Figure  5-  By  differentiating  (12)  with 
respect  to  r and  using  Eq.  (6),  wo  obtain 


r1d(y1  - y12)2  " r2d(y2  ~ y12)2 


(KlL  r.  dr2 


(13) 


Substituting  (7)  and  (8)  yields  the  simple  result 


U-d<y1  - y12)2  * U(r)dr2 


(14) 
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When  y2  is  set  equal  to  b/2,  Eq.  (1*!)  is  valid  for  a tip 
roll-up.  When  U(r)  « , the  results  given  in  Eqs.  (4)  and  (9) 

are  readily  obtained 


Befor®  coupling  the  wing  drag  distribution  to  the  axial 
velocity  In  the  vortex,  it  is  possible  to  show  how  the  axial 
velocity  modifies  the  velocity  V(0)  at  the  center  of  the  vortex, 
Assuming  that  U(0)  is  finite,  Eq . (14)  may  be  integrated  for 
small  r to  yield 

r n "1  i/2 


" u.  " 

uToTj 
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(y12  - yx) 


as  v + 0 . Following  Ref.  13,  the  tangential  velocity  for  the 
"interior*  vortex  is 

r * 1 (y±)  “ ^ ^ p ^ 

■utr)  ^ 1— L£ l * ± £—  (- 

V NX  > On-™  V * 


As  r •+•  0 , we  may  write 

V(0)  ,hr£\ 
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ayj\ 

. 1 dr 
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Calculating  av^'dr  and  dy2/dr  from  Eqs.  (15)  ftnd  (8)  and 
substituting  these  values  into  Eq.  (17)  yields 


Y<0) 
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For  a tip  vortex,  Eq.  (15)  becomes 

r u i1/2 

QO 

r “ uToT 

and  (16)  is  written 
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(19), 

taking 

the 

yields 

dr  b/2 
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For  both  "interior”  ar.d  tip  vortices,  deficits  in  axial  velocity 
(U(0)/U-  < 1)  result  in  a reduction  of  the  inviscid  centerline 

swirl  velocity. 
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We  may  now  turn  our  attention  to  coupling  the  axial  velocity 
in  the  vortex  to  the  wing  drag  distribution.  Making  an  axial 
momentum  balance  across  a "cylindrical*’  control  volume  of  radius 
r containing  the  portion  of  the  wing  between  y^  and  y„,  yields 

C y2  f* 

I d(n)dn  L 2it  [p  + pU(U  - Uw)]cd5  = 0 (22) 

J Jo 

where  the  axial  velocity  of  the  fluid  fluxing  through  the  "cylin- 
drical" surface  is  approximated  by  U . The  pressure  far  upstream 
has  been  taken  to  be  zero  and  d is  the  wing  sectional  drag.  When 
the  U2  term  is  linearized,  Eq.  (22)  is  that  given  by  Brown  (Ref. 
10).  Equation  (22)  is  written  in  the  same  spirit  as  Eqs . (1)  - (5) 
since  it  assumes  that  the  wing  drag  distributes  itself  in  the 
rolled-up  vortex  in  the  same  manner  as  the  shed  axial  vorticity. 

As  discussed  by  Brown,  the  assumption  is  a natural  one  since  the 
axial  vortex  lines  and  the  viscous  wake  are  one  and  the  same. 

Differentiating  Eq.  (22)  and  substituting  (14)  yields 

dy?  ["  1 d(y.  - y1?)2 

d(yx)  - d(y2)  ♦ p(U  - U.)J  ^ — (23) 

For  a tip  roll-up,  dy2/dy^  ■ 0.  The  nonuniform  pressure  in  the 
vortex  is  primarily  a result  of  the  swirl  and  may,  therefore,  be 
calculated  from 


P * - I ^ dc  (24) 

4*2  Jr  t3 

Equations  (14),  (23),  and  (24)  form  the  nonlinear  system  to  be 
solved.  The  boundary  conditions  are 


P 


(25) 


Together,  Eqs.  (14),  (23),  (24)  and  (25)  determine  the  inviscid 
axial  and  swirl  velocities  in  either  an  "interior"  or  tip  vortex 
given  the  lift  and  drag  distribution  over  the  appropriate  region 
of  the  wing  from  which  the  vortex  develops.  Nonlinearity  and  the 
nature  of  the  boundary  conditions  dictate  that,  in  general,  solu- 
tions will  have  to  be  obtained  numerically. 
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At  this  point,  it  is  appropriate  to  find  a simple  analytical 
solution  which  will  illustrate  the  effect  of  wing  drag  distribu- 
tion on  inviscid  vortex  structure.  Unfortunately,  the  direct 
problem  of  specifying  the  lift  and  drag  distribution  does  not 
appear  to  yield  analytic  solutions.  However,  the  indirect  problem 
of  specifying  the  axial  velocity  and  lift  distribution  and  deter- 
mining the  drag  distribution  and  swirling  velocity  is  straightfor- 
ward for  a simple  distribution. 


Equation  (14)  can  be  integrated  if  the  wing  loading  is 
linear  (a  tip  roll-up)  and  therefore  of  the  form 


r - r6(i  - £)  (26) 

and  the  axial  velocity  is  given  by 


i * 

2*1 

U « U„  a + b ( j 

=?)] 

(27) 

The  constants  a and  b may  be  chosen  such  that  U is  positive; 
therefore,  a > 0.  Negative  axial  velocities  imply  a flux  of 
angular  momentum  from  downstream  and  violate  the  assumptions 
implied  in  Eq.  (12).  The  radius  of  the  region  containing  all  the 
vorticity  shed  between  the  wing  root  and  tip  rt  is  to  be  deter- 
mined. Integrating  Eq.  (14)  yields 


(s 


(28) 


When  r ■ rfc  , y 


0 


and,  therefore, 


(29) 


The  Betz  result  is  obtained  with  a * 1,  b = 0;  all  the  vorticity 
is  contained  within  a circle  having  radius  b/4.  Referring  to 
Eq.  (27),  sufficient  conditions  for  an  axial  velocity  excess  in 
the  wake  occur  when  a > 1 and  b >0;  rt  decreases  and  the 
vortex  is  intensified  in  that  the  swirl  velocities  are  increased. 
Axial  velocity  defects  are  associated  with  increases  in  r^  and, 
therefore,  more  diffuse  vortices.  It  is  important  to  remember, 
however,  that  since  the  lift  distribution  is  unchanged  the  total 
axial  flux  of  angular  momentum  from  the  region  containing  axial 
vorticity  is  not  changed.  Therefore,  the  intensity  of  the  vortex, 
as  measured  by  the  flux  of  angular  momentum,  is  unchanged  by  drag. 
The  deintensification  which  does,  in  fact,  occur  is  brought 
about  by  redistributing  the  angular  momentum  outward  so  that  small 
encountering  aircraft  could  Interact  with  less  of  the  vortex. 

Since  drag  is  likely  to  result  in  higher  turbulence  levels  in  the 
vortex  (at  least  initially),  it  is  likely  that  the  outward 
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redistribution  of  angular  momentum,  as  calculated  here,  is  further, 
aided  by  turbulent  processes. 


The  circulation  and  swirl  velocity  distributions  in  the 
vortex  are  given  by 


(30) 


(31) 


Referring  to  Figure  6,  the  drag  on  the  wing  is  calculated  from 
Eq.  (23)  for  five  cases.  At  an  average  drag  coefficient  of  about 
0.012  (a  typical  value),  the  axial  velocity  in  the  vortex  is 
uniform  and  equal  to  the  free  stream  value  (case  3).  The  vortex 
radius  is  0.55  b/2  and  is  taken  to  define  a reference  circular 
area  so  that  we  may  illustrate  the  deintensification  which  results 
as  a consequence  of  wing  drag.  By  increasing  the  average  drag 
coefficient  by  nearly  an  order  of  magnitude  to  0.11,  the  vortex 
radius  is  increased  16J  to  0.585  b/2  (case  2).  The  flux  of  angu- 
lar momentum  through  the  reference  area  when  compared  with  case  3 
is  one  measure  of  the  deintensification  which  can  be  achieved  by 
increased  drag.  The  calculation  shows  that  the  flux  In  case  2 is 
reduced  a highly  desirable  43*.  However,  the  model  suggests  that 
vortex  deintensification  brought  about  solely  by  this  technique 
will  have  prohibitively  large  drag  penalties.  The  problem  of  vhe 
least  intense  vortex  for  a given  lift  and  drag  coefficient  is 
surely  worthy  of  additional  study. 


It  is  now  apparent  that  axial  velocities  in  the  vortex  may 
be  either  an  excess  or  defect  depending  on  the  wing  lift  and  drag 
distribution.  In  addition,  it  is  quite  possible  that  axia.l 
velocity  distributions  can  result  which  have  an  excess  over 
certain  radial  intervals  and  defects  over  others.  Figure  7 
illustrates  such  a situation.  Tailoring  the  drag  distribution  so 
that  strong  axial  velocity  gradients  develop  may  prove  to  be  an 
effective  means  of  enhancing  turbulent  decay. 


To  summarize  the  results  of  this  section,  it  has  been  shown 
that  the  Bets  model  and  extensions  to  calculate  the  roll-up  of 
aircraft  vortices  may  be  further  modified  to  include  the  effect 
of  wing  drag  distribution  on  vortex  structure.  Simple  sample 
calculations  show  that  increased  drag  reduces  the  intensity  of 
the  vortex  by  redistributing  the  flux  of  angular  momentum  over  a ' 
greater  area  in  the  wake.  Large  restructuring  of  the  wake  by 
increased  drag,  however,  appears  to  require  large  drag  penalties. 


3.  ROLL-UP  OF  THE  TURBULENT  BOUNDARY  LAYER 

Having  developed  a model  to  predict  the  mean  swirl  and  axial 
velocity  distributions  in  wake  vortices,  there  remains  the  problem 
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of  estimating  how  the  turbulent  kinetic  energy  from  the  wing 
distributes  itself  in  the  rolled-up  vortex.  Once  this  is  known* 
it  is  possible  to  make  detailed  calculations  of  the  local  turbu- 
lent decay  of  these  vortices  using  a computer  program  which  was 
developed  by  A.R.A.P.  (Ref.  15)  for  ARL-  One  might  then  ask  the 
questions  as  to  how  much  turbulence  is  required  and  where  on  the 
wing  would  we  introduce  it  to  bring  about  the  most  rapid  aging  of 
the  vortex.  The  answers  to  these  questions  have  obvious  applica- 
tions. We  will  only  treat  the  roll-up  of  a tip  vortex  here, 
although  extending  the  analysis  to  interior  vortices  is  straight- 
forward . 

In  order  to  derive  a simple  expression  which  gives  the 
distribution  of  turbulent  kinetic  energy  q*z(r)  in  the  rolled-up 
vortex,  it  is  necessary  to  make  a somewhat  tenuous  assumption; 
namely,  the  production,  dissipation,  and  diffusion  of  turbulence 
are  negligible  during  the  roll-up  process.  We  might  refer  to  the 
turbulence  as  being  in  a "frozen”  state  so  that  we  may  write 

C +h(y )/2  P ? 

I q*(y,z)U(y ,z)dz  dy  * - q'^(r)2irrU(r)dr  (32) 

J-h(y)/2 


Equation  (32)  states  that  the  flux  of  turbulent  kinetic  energy 
immediately  behind  the  wing  through  the  area  h(y)dy  equals  the 
flux  in  the  annular  area  2nr  dr.  h(y)  must  be  calculated  from 
mass  conservation.  Substituting  Eq.  (14)  into  (32)  yields 


q,2(r) 


-P(y>/V< y,z>  M 

J-  h(y)/2 


U_ 


- ^(y) 


(33) 


— 2 

h.Cy)  * itd(y  - y)  /dy  and  is  the  height  of  a rectangular  area  of 
width  dy  through  which  the  volume  flux  equals  that  in  the 
annular  area  2trr  dr.  The  relationship  between  h.  and  h is 
simply  1 


h1(y) 


/. 


+h(y)/2 


-h(y )/2 


^^-dz 


(31*) 


• and  h accounts  for  the  mass  defect  in  the  turbulent  boundary 

layer.  Since  h is,  in  general,  much  greater  than  the  turbulent 
[ boundary  layer  thickness  (except  at  the  wing  tip),  the  integral 

| in  (33)  is  independent  of  h(y)  and,  as  the  vortex  rolls  up,  the 

? turbulent  kinetic  energy  in  the  vortex  will  diminish  as  h.(y)*'-1- 

l The  drop-off  outward  from  the  vortex  center  is  a consequence  of 

the  fluid  containing  turbulent  kinetic  energy  near  the  wing  root 
",  being  distributed  over  a larger  area  in  the  vortex  than  that 
' fluid  near  the  wing  tip. 
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4.  INTERACTIVE  MOTION  OP  SHED  VORTICITY 


To  complete  the  description  of  an  aircraft  wake,  there 
remains  the  problem  of  determining  how  the  discrete  vortices, 
whose  structure  has  been  determined  with  the  models  developed  in 
the  previous  sections,  move  under  each  other’s  influence.  We  are 
not  concerned  here,  however,  with  self-induced  motions  such  as 
those  which  in  the  presence  of  other  vortices  give  rise  to  the 
urow  instability  (Ref.  16).  Our  aim  here  is  only  to  explore  wake 
geometry  by  determining  the  approximate  locations  of  the  discrete 
vortices  which  form  the  wake. 

We  start  with  what  is  now  the  commonplace  practice  of  model- 
ing discrete  vortices  with  infinite  two-dimensional  irrotational 
vortices.  This  assumption  is  correct  as  long  as  gradients  in  the 
flight  direction  of  the  wake  are  small.  When  this  condition  is 
met,  the  radii  of  curvature  of  the  vortex  filaments  are  large  and 
the  self-induced  motion  may  be  neglected.  The  motion  of  the 
vortices  then  does  not  depend  crucially  on  the  vortex  structure 
(see  Ref.  11).  The  condition  for  small  axial  gradients  is  one 
which  requires  that  CL/AR  be  small,  since  the  motions  perpendi- 
cular to  the  iiight  direction  proceed  with  velocity  of  order  f/b 
and  axial  gradients  are,  therefore,  of  order  r/bU®  . Obviously 
this  condition  is  met  under  most  circumstances  since  - 6^(1) 
and  AR  is  typically  0’(1O).  We  may  then  make  the  equivalence 
between  time  and  downstream  distance  through  d/dt  = Uoo(d/dx). 

We  consider  the  motion  of  n irrotational  vortices  of 
strength  = T±/2t\  at  instantaneous  locations  (y^,  f’-e 

velocity  of'  motion  of  tne  jth  vortex  is  equal  to  the  fluid  motion 
at  the  point  (y«,  z*)  in  the  absence  of  the  Jth  vortex.  Mathe- 
matically, the  instantaneous  velocity  is  calculated  from 


(35) 


(36) 


and  the  summation  is  to  exclude  the  jth  vortex.  r^..  is  the  dis- 
tance between  the  jth  and  ith  vortex  and  is  given  By 

ri3  ■ (i,i  - »3)2  + <*i  - zj)2 

Figure  8 illustrates  the  configuration  to  be  studied.  Some  sim- 
plification now  results  since  the  vortices  occur  in  pairs,  and  it 
is  not  difficult  to  show  from  the  calculation  of  the  vertical 
momentum  of  the  flow  that 
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where  J is  a constant.  In  general,  however,  except  for  a few 
special  cases,  the  solution  to  Eqs.  (35)  and  (36)  must  be  obtained 
by  numerical  integration. 


If,  for  the  time  being,  we  do  not  concern  outselves  with  the 
time  history  of  the  motion,  we  may  then  obtain  information 
regarding  the  relative  trajectories  of  the  vortices.  Dividing 
Eq.  (35)  by  (36),  it  is  quite  remarkable  that  the  expression  can 
be  integrated  to  yield 

n n 

w » - l l k,k.  log  r.,  (39) 

1 j 1 J 
(i*J) 


It  is  readily  recognized  that  W is  the  integral  by  noting  that 
(35)  and  (36)  may  be  written  in  the  form 
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dyj  _ 3W_ 
dt  " 3z^ 
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(40) 

(41) 


It  may  be  shown  that  W Is  related  to  that  portion  of  the  kinetic 
energy  of  the  flow  that  depends  on  only  the  relative  distance 
between  vortices. 


While  the  aircraft  wake  may  be  made  up  of  any  number  of  pairs 
of  vortices,  it  i3  likely  that  nu_st  aircraft  wakes  can  be  described 
by  three  or  fewer  pairs.  When  two  pairs  are  adequate,  a great  deal 
can  be  said  regarding  the  geometry  of  the  wake  without  resorting 
to  numerical  calculation.  For  two-pair  wakes  having  a common  axis 
of  symmetry  (the  aircraft  centerline),  Eq.  (39)  can  be  written  in 
the  form 


B 


Le2  + (y,  - y3rJ 


(C./k,  K-/k 

) 3 x(y.)  1 


(42) 


where  e is  1/2  the  vertical  separation  distance  between  pairs. 
yl  or  y3  may  be  eliminated  from  (42)  by  substituting  (38).  The 
constant  B is  to  be  evaluated  with  e * 0 (the  trailing  vortex 
pairs  lie  In  the  same  plane  at  the  wing).  Equation  (42)  then 
describes  all  possible  relative  vortex  trajectories  with  a 
rather  simple  expression.  The  obvious  question  is,  "For  what 
values  of  J and  B are  trajectories  which  have  e -*■  °» 
possible?"  The  investigation  of  the  relative  trajectories,  while 
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very  interesting,  is  tedious  and  we  only  present  the  results  here. 

It  is  possible  to  construct  what  might  be  called  a wake  classifi- 
cation chart  for  four  vortex  wakes,  as  shown  in  Figure  9*  Prom 
this  chart,  knowing  the  relative  vortex  strengths  and  Initial 
locations.  It  is  possible  to  determine  If  the  wake  geometry  will 
be  such  that  the  vortex  pairs  will  remain  contiguous  or  separate 
downstream. 

We  have  carried  out  four  sample  calculations,  shown  in 
Figures  10-13.  The  trajectories  are  viewed  from  downstream  behind 
the  aircraft.  The  quite  diverse  behavior  of  the  vortices  is 
predicted  by  Figure  9 where  we  have  shown  where  '.hese  cases  lie  on 
the  wake  classification  chart.  The  ability  to  predict  wake  behavior 
should  prove  to  be  of  great  use  to  the  aircraft  designer  as  well  as 
to  the  aerodynamicist  concerned  with  wake  alleviation. 

When  the  aircraft  vortex  wake  must  be  modeled  by  three  or  more 
pairs  of  vortices,  trajectories  can  only  be  obtained  by  direct 
numerical  calculation.  To  check  our  numerical  scheme  for  three- 
pair  wakes,  we  chose  to  model  the  wake  of  an  aircraft  where  flow 
visualization  was  available,  larlett  and  Shivers  (Ref.  18) 
describe  such  a test  of  an  E.B.F.V/STOL  model.  The  aircraft  was  one 
which  had  a blown  inboard  flap  and  was  consequently  highly  loaded 
there.  From  data  presented,  we  crudely  estimated  the  lift  distri- 
bution due  to  circulation  lift  to  be  that  shown  in  Figure  14.  In 
Figures  15-17,  the  results  of  the  numerical  solution  of  Eqs . (35) 
and  (36)  are  presented.  The  aircraft  is  at  10°  angle  of  attack 
and  Is  a scale  representation  of  the  test  model.  The  results  are 
shown  In  perspective  to  help  resolve  the  complicated  wake  geometry 
that  results.  As  can  be  seen,  a strong  outboard  flap  vortex-tail 
interaction  is  predicted  by  the  calculation.  Figure  18  has  been 
reproduced  from  Figure  17b  in  Reference  18  and  clearly  shows  the 
interaction  Just  described.  It  is  reported  in  Reference  18  that 
the  aircraft  tested  became  longitudinally  unstable  at  angles  of 
attack  in  excess  of  7°  due  to  the  strong  downwash  field  induced  at 
the  tail.  Simple  calculations  of  this  type  should  prove  to  be 
extremely  useful  In  that  they  are  able  to  predict  the  observed 
unfavorable  vortex  tall  interference.  With  an  estimate  of  the 
rate  at  which  vortex  sheets  roll  up  (a  simple  model  is  developed 
to  do  this  In  Appendix  A)  and  with  the  models  previously  described, 
the  ability  to  accurately  predict  the  downwash  field  In  the  near 
wake  will  be  possible. 
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classification  chart  for  two-vortex-pair  wakes 


TRAILER  TRAJECTORIES 
OUTBOARD  FLAP 
sf/s=.33  8»=  40°  r,/r,  =-l.l 
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Figure  12.  Vortex  centroid  locations  as  seen  from  downstream 
pairs  of  opposite  sign.  Pairs  remain  together. 
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TRAILER  TRAJECTORIES 
INBOARD  FLAP 

Sf/s  =.65  8f  = 40° 

rt/rf  =-.35 


Vortex  centroid  locations  as  seen  from  downstream; 
pairs  of  opposite  sign.  Weak  interaction  between 
pairs;  pair3  divergrt. 


Figure  15.  Calculated  wake  geometry;  top  view 
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gure  16.  Calculated  wake  geometry;  aide  elevation 


OBLIQUE  VIEW 


Figure  17.  Calculated  wake  geometry;  ob'.ique  view 


Figure  18.  Smoke  visualization  showing  the  outboard  flap  vortex 
in  proximity  to  the  stabilizer 
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SECTION  III 

THE  EXPERIMENTAL  EFFORT 


The  validity  of  the  Betz  models  and  those  extensions  developed 
in  this  report  were  checked,  where  possible,  with  experimental 
measurements  made  in  the  V/STOL  tunnel  at  NASA  Langley  Research 
Center.  This  section  describes  the  design  and  implementation  of 
the  test  program  which  was  carried  out. 


i 


1.  WAKE  CLASSIFICATION  AND  EXPERIMENT  DESIGN  . 

In  order  to  check  the  models  which  describe  the  mean  wake 
structure  (th^  turbulent  structure  was  not  to  be  measured  here),  it 
was  necessary  to  obtain  wing  lift  and  drag  distributions  and  the 
velocity  distributions  in  the  downstream  wake.  The  test  model  and 
tunnel  facility  were  supplied  by  the  government. 

It  was  decided  that  the  test  would  concern  Itself  with  two- 
vortex-pair  wakes.  One  pair  would  result  from  an  "interior"  roll- 
up and  the  second,  obviously,  as  a result  of  a tip  roll-up.  The 
results  of  the  two-vortex-pair  trajectory  model  suggest  that  the 
geometry  of  two-vortex-pair  wakes  may  be  classified  such  that  they 
fall  Into  one  of  four  categories: 

1)  Both  pairs  are  of  the  same  sign  and  the  pairs  remain 
together  (e.g.,  see  Fig.  10). 

2)  Both  pairs  are  of  the  same  sign  and  the  pairs  separate 
(e.g. , see  Fig.  11) . 

3)  The  pairs  are  of  opposite  sign  and  remain  together 
(e.g. , see  Fig.  12) . 

4)  The  pairs  are  of  opposite  sign  and  separate 
(e.g.  , see  Fig  13) . 

Categories  (3)  and  (4)  were  Judged  unlikely  to  occur  under  current 
aircraft  operating  conditions  and,  therefore,  only  one  configura- 
tion which  produced  two  vortex  pairs  of  opposite  sign  was  tested. 

The  model  supplied  by  the  government  is  shown  with  relevant 
dimensions  in  Figure  19.  Six,  roughly  equal  plain  flaps  span  the 
entire  trailing  edge  of  the  wing.  Flap  deflection  angles  were 
fixed  by  brackets  but  covered  a wide  range  so  that  the  number  of 
configurations  possible  was  very  large.  With  one  configuration 
being  the  clean  wing,  the  four  test  configurations  are  shown 
schematically  in  Figure  20. 

The  tunnel  facility  was  the  V/STOL  tunnel  at  NASA  Langley. 

The  tunnel  test  section  measures  14.5  ft  high,  21.75  ft  wide,  and 
50  ft  long.  The  tunnel  speed  range  is  0-200  knots.  The  model  was 
sting-mounted  as  shown  in  Figure  21.  During  testing,  the  model 
was  positioned  as  near  to  the  tunnel  centerline  as  possible. 
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Figure  20.  Model  configurations  tested  in  NASA  Langley  V/STOL 
tunnel 
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Figure  21.  Test  model 
sting-mounted  in  the 
V/STOL  tunnel  at  NASA 
Langley  Research  Center 


Figure  22.  Model  test 
configuration  4 


Figure  23.  The  wake 
traversing  mechanism 
mounted  with  the  hot- 
fllm  probe  5 chords 
downstream  of  the 
trailing  edge 


Downstream  wake  velocities  were  measured  using  a Thermal  Systems 
Model  1050,  split-film,  three-component  anemometer  mounted  on  a 
traversing  mechanism  shown  in  Figure  23.  Tunnel  test  speed  was 
85  ft'e-ec  (dynamic  pressure  of  9 lb/ft* ) to  avoid  structural 
problems  with  the  traversing  mechanism. 


Wing  load  and  drag  distributions  are  the  inputs  to  the  theo- 
retical models  to  determine  wake  structure.  Load  distribution  was 
determined  from  integration  of  pressure  data  from  existing  taps 
located  at  six  spanwise  locations.  Because  of  the  importance  of 
these  data,  NASA  Langley  personnel  also  performed  lifting  surface 
calculations  on  the  geometries  to  be  tested  to  provide  a check  on 
the  measured  distributions.  Profile  drag  distribution  data  were 
obtained  by  measuring  the  axial  velocity  defect  at  the  wing 
trailing  edge  and  determining  the  drag  from  an  axial  momentum 
balance.  Since  lift  distribution  data  would  only  be  given  by  six 
spanwise  data  points,  vertical  traverses  to  determine  the  axial 
velocity  defect  were  only  made  at  the  same  six  spanwise  stations. 

To  assure  a turbulent  boundary  layer  on  the  upper  side  of  the  wing, 
a trip  strip  (0.06-inch  diameter  particles)  0.1-inch  wide  was 
placed  along  the  wing  1 inch  from  the  leading  edge. 

Prior  to  fixing  the  test  model  angle  of  attack  and  making 
velocity  measurements  in  the  wake,  flow  visualization  studies  were 
undertaken  to  determine  if  the  chosen  model  configurations  pro- 
duced the  desired  wake  structure.  The  trailing  vortices  were 
marked  with  smoke  introduced  upstream  of  the  wing.  The  smoke, 
upon  passing  around  the  wing,  becomes  rolled  up  with  the  fluid 
forming  the  trailing  vortices.  Results  from  prior  NASA  Langley 
test  programs  have  shown  that  this  technique  gives  a very  good 
description  of  the  vortex  positions. 
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Velocity  measurements  in  the  wake  were  to  be  made  at  the  wing  j 

trailing  edge,  5 chords,  and  30  chords  downs  am.  However,  the 
30-chord  measurement  was  subject  to  being  mov  u closer  to  the  wing 
if  vortex  meander  was  Judged  so  severe  that  meaningful  data  could 
not  be  obtained.  V/STOL  tunnel  personnel  have  had  considerable  , 

experience  in  making  these  measurements.  Because  of  the  basic  \ 

unsteadiness  in  the  wake,  the  anemometer  signal  was  sampled  over  a j 

t time  interval  of  5 seconds  and  stored  on  magnetic  tape.  Prior  test 

| experience  has  shown  that  this  sample  time  interval  is  a good  t 

I compromise  between  obtaining  meaningful  averages  and  using  excess- 

i ive  tape  storage  capacity.  A disadvantage  of  this  system  is  that 

no  on-line  capability  was  available,  and  first  examination  of  the 
velocity  data  could  not  be  made  until  after  testing  was  complete. 

The  procedure  used  to  determine  the  approximate  location  of 
the  vortex  center  downstream  was  to  mount  an  impact  probe  above  the 
hot-film  sensor  and  scan  until  a minimum  in  total  pressure  was 
; obtained.  Since  the  vortex  is  composed  of  vortical  material  shed 

; from  the  boundary  layer  of  the  wing,  the  head  in  the  vortex  must  be 

• less  than  the  surrounding  irrotational  flow.  To  aid  in  determining 

* the  area  to  be  scanned,  smoke  was  again  introduced  upstream  of  the 


\ 
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wing  which  rolled  up  and  marked  the  vortex  center.  The  hot-film 
sensor  was  shielded  to  avoid  contamination  and  the  impact  probe  was 
visually  lined  up  with  the  vortex  center. 

At  the  5-chord  and  30-chord  downstream  positions,  three  hori- 
zontal traverses  were  made  across  each  vortex.  The  traverses  were 
spaced  1 inch  vertically  apart,  and  velocity  data  were  obtained  at 
approximately  0.5-inch  horizontal  intervals  across  the  vortex. 

This  grid  spacing  provides  a good  description  of  the  vortex  struc- 
ture while  not  requiring  a prohibitively  large  number  of  data 
points.  Traverses  at  tne  wing  trailing  edge  were  made  vertically, 
the  distance  between  data  points  being  approximately  0.1  inch. 


2.  THE  TEST  SEQUENCE 

The  test  program  was  performed  in  the  following  sequence. 

1)  Force,  moment,  and  wing  pressure  data  were  obtained  for 
the  four  model  configurations  through  an  angle  of  attack 
range  which  would  bracket  the  estimated  test  angle  of 
attack.  Pressure  data  were  oDtained  prior  to  fixing  the 
angle  of  attack  after  flow  visualization  studies  to  avoid 
plugging  the  small  pressure  orifices  with  oil  residue 
from  the  smoke  generators . 

2)  Flow  visualization  was  carried  out  using  smoke  and  helium- 
filled  soap  bubbles  to  mark  the  wake  flow  field.  Photo- 
graphs recorded  the  trailer  positions  from  cameras  which 
were  permanently  mounted  and  hand  held.  Test  angles  of 
attack  were  fixed  upon  obtaining  the  desired  wake 
geometries. 

3)  Velocity  measurements  were  made  in  the  wake  using  the 
hot-film  sensor  mounted  on  an  existing  traversing  mechan- 
ism which  was  controlled  outside  the  tunnel  test  section. 
Probe  calibration  was  checked  prior  to  and  at  the 
completion  of  each  test  run. 


SECTION  IV 


TEST  RESULTS  AND  COMPARISON  WITH  THEORETICAL  MODELS 


The  model  test  configurations  and  test  conditions  are 
summarized  below. 


Configuration 

Flaps 

a0  Re„ 

___  c_ 

q,  lb/ft 

1 

None 

- 

10  2.25  x 105 

9 

2 

Inboard 

30 

6 j 

3 

Inboard  & 

Midspan 

15 

6.5  | 

4 

Outboard 

30 

1.5  ^ 

' 

’ 

The  results  of  the  test 

program 

are 

1)  Wing  lift  distribution  obtained  from  integration  of 
pressure  tap  measurements  for  the  four  model  test 
configuration; ; 

2)  Photographs  of  the  wake  structure  for  the  four  model 
test  configurations; 

3)  Three-component  mean  velocity  measurements  in  the  wake 
at  the  wing  trailing  edge  and  5 chords  downstream  for 
the  four  model  test  configurations.  Three-component 
mean  velocity  measurements  in  the  wake  10  chords  down- 
stream for  two  model  test  configurations.  Wind  tunnei 
scheduling  problems  prohibited  testing  two  additional 
configurations  at  the  10-chord  position. 

1.  WING  PRESSURE  DATA  AND  WING  LIFT  DISTRIBUTIONS 

Wing  pressure  data  were  supplied  by  FASA  Langley  Research 
Center  in  pressure  coefficient  form.  Wing  lift  distribution  was 
obtained  by  integrating  the  pressure  data  over  the  chord  of  the 
airfoil  and  then  taking  the  vertical  component  of  the  resulting 
force  coefficient.  In  general,  pressure  data  on  deflected  flaps 
were  difficult  to  interpret  due  to  the  limited  number  of  pressure 
taps  (two)  on  the  underside  of  the  flap.  When  this  problem 
occurred,  data  on  the  flap  were  faired  to  yield  what  were  judged 
to  be  reasonable  chordwise  distributions. 

The  experimentally  determined  wing  lift  distributions  are 
shown  in  Figure  24.  The  limited  number  of  spanwise  pressure  tap 
locations  forced  some  fairing  of  the  results.  The  distributions 
are  physically  very  reasonable,  however,  and  will  produce  the 
desired  wake  structures.  All  configurations  had  a significant- 
lift  drop-off  inboard  as  recorded  by  the  pressure  taps  at 
2y/b  m.  0.15.  This  is  a result  of  the  rather  large  model  fuselage 
which  is  needed  to  hold  instrumentation.  The  strongly  loaded  tip 
in  the  clean  configuration  is  believed  to  arise  as  a consequence 
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Figure  24,  Lift  distributions  from  wing  surface  pressure 
measurements 
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of  vortex  lift  much  like  that  which  results  on  a delta  wing. 

Smoke  flow  visualization  of  this  configuration  indicates  that  the 
flow  separates  at  the  square  tip,  and  roll-up  actually  proceeds 
from  the  wing  tip  near  the  leading  edge. 

The  experimentally  determined  lift  distributions  may  be 
compared  with  lifting  surface  calculations  supplied  by  the  V/STOL 
tunnel  staff.  The  calculations  are  shown  in  Figure  25  and  match 
model  geometry  and  measured  lift  coefficient.  It  is  obvious  that 
there  are  some  differences  between  these  results.  Both  results 
are  subject  to  error,  and  it  is  likely  that  the  true  lift  distri- 
bution Is  somewhere  between  the  two  results.  We  view  both  results 
cautiously  for  several  reasons.  On  one  hand,  the  lifting  surface 
calculation  neglected  the  rather  large  fuselage.  Perhaps  more 
jr  significant  is  the  fact  that  the  test  Reynolds  number  was  moderate 

and  the  airfoil  section  had  a thickness  ratio  of  0.18.  On  the 
other  hand,  however,  experimentally  determined  distributions,  as 
previously  mentioned,  suffered  from  a limited  number  of  pressure 
taps.  So  that  no  bias  would  be  shown  in  the  comparison  of  measured 
wake  velocity  fields  with  those  calculated  from  wing  lift  and  drag 
distributions,  calculations  were  performed  using  both  experimental 
and  lifting  surface  results. 


2.  FLOW  VISUALIZATION 

Results  of  the  flow  visualization  studies  ari  calculated  vortex 
positions  are  shown  In  Figures  26-33.  The  vortices  were  made 
visible  with  smoke,  as  previously  described.  Vortex  centroid  "osi- 
tions  were  calculated  using  the  model  developed  in  Section  II. 4 and 
are  shown  on  consecutive  figures.  As  an  Input  to  the  calculation, 
the  lifting  surface  load  distributions  were  divided  into  flap  and 
tip  vortices  according  to  the  criterion  previously  described.  For 
the  small  downstream  distances  observed  here,  the  results  are  quite 
insensitive  and  differ  little  from  the  results  obtained  using  the 
experimentally  determined  lift  distributions.  The  calculated 
trajectories  are  presented  in  perspective  with  the  viewing  position 
taken  to  be  the  camera  location.  As  can  be  seen,  the  agreement  Is, 
In  general,  very  good. 

Photographs  taken  with  cameras  permanently  mounted  in  the 
tunnel  provided  somewhat  less  information  than  those  taken  with 
hand  held  cameras.  In  Figure  3^  are  shown  side  views  of  the  wakes 
of  configurations  2 and  4 taken  with  a permanently  mounted  camera. 

Helium  filled  soap  bubbles  were  also  introduced  upstream  of 
the  wing  In  an  effort  to  visualize  the  downstream  wake.  Although 
the  bubbfcs  did  roll  up  into  the  vortex,  observation  was  hampered  by 
the  low  density  of  bubbles  and  difficulties  with  illumination. 

3.  COMPARISON  OF  COMPUTED  AND  MEASURED  WAKE  VELOCITY  DISTRIBUTIONS 

Wake  velocities  were  received  in  Cartesian  tunnel  coordinates, 
positive  x measured  downstream  from  the  trailing  edge  and  positive 
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Figure  28.  Smoke  visualization;  configuration  2.  (a)  viewed  from 

downstream;  (b)  viewed  from  downstream  and  above 
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y measured  from  the  tunnel  centerline  to  the  right  wing  tip 
(facing  upstream).  Positive  z is  measured  vertically  upward 
from  the  tunnel  floor.  Velocity  measurements  were  made  behind  the 
left  half  of  the  wing,  the  wing  tip  being  located  at  y « -40  in. 

The  /elocity  traverses  made  at  the  wing  trailing  edge  (actually 
within  3 inches  of  the  trailing  edge)  were  reduced  to  estimate  the 
profile  drag  of  the  wing.  Profile  drag  distributions  were  calcula- 
ted from  the  axial  velocity  deficits  shown  in  Figures  35  and  36 
for  configurations  1 and  2,  respectively.  The  computed  sectional 
drag  coefficients  are  shown  in  Figure  37.  The  large  value  of 
profile  drag  at  the  inboard  stations  for  the  clean  configuration 
suggests  that  the  flow  is  beginning  to  separate. 

NASA  tunnel  staff  members  have  confirmed  (with  tuft  studies) 
separated  flow  at  the  unfaired  fuselage  wing  junction  for  this 
model  at  moderate  angles  of  attack. 

Wake  velocity  measurements  at  distances  further  than  10  chords 
downstream  were  precluded  by  severe  wake  meander.  Vortex  centers 
marked  by  smoke  had  excursions  whose  amplitudes  were  visually 
estimated  to  be  nearly  1 foot  at  downstream  distances  of  30  chords. 
The  validity  of  averaged  velocity  data  under  such  circumstances  wa3 
doubted,  and  measurements  were  made  at  5 and  10  chords  downstream 
where  meander  was  judged  to  be  moderate.  Experimental  velocity 
data  are  shown  in  Figures  38-61  and  Figures  B.1-B.4.  At  each  down- 
stream station,  three  horizontal  traverses  were  made  at  three 
elevations  across  each  vortex.  The  traverses,  in  principle,  were 
to  pass  above,  through,  and  below  the  vortex  center,  a task  not 
easily  achieved. 

Since  three  traverses  were  made  across  each  vortex,  it  was 
decided  to  check  the  Betz  model  and  extensions  by  trying  to  dupli- 
cate the  velocities  measured  with  each  traverse  by  computation. 
Comparison  in  this  way  is  a more  stringent  check  of  the  models  than 
has  previously  been  used  in  the  literature  where  it  is  often  the 
practice  to  U3e  the  velocity  measurements  from  orly  one  traverse 
and  to  fit  tne  computed  results  to  the  experimental  data  oy  adjus- 
ting the  location  of  the  vortex  center.  Even  small  changes  in  the 
location  of  the  vortex  center  can  result  in  large  changes  in  the 
velocities  computed  at  a point  due  to  the  large  gradients  which 
occur  near  the  vortex  center.  When  two  or  more  traverses  at  differ- 
ent elevations  are  available,  the  freedom  to  choose  vortex  center 
locations  is  limited  in  that  only  one  location  will  distribute  the 
error  in  a prescribed  manner. 

Swirling  velocity  distributions  in  the  wake  were  calculated 
from  the  Betz  model  for  tip  vortices  and  with  the  extension  by 
Ref.  13  for  flap  vortices  (wing  drag  is  neglected  here).  The  wing 
load  distributions  which  would  result  in  multiple  vortex  wakes  were 
divided  according  to  the  criterion  given  in  Ref,  13  and  in  Section 
II. 1.  A computer  program  was  then  developed  to  best  fit  the 
experimental  velocity  measurements  with  the  computed  values.  The 
procedure  was  to  minimize  the  square  of  the  difference  between  the 
experimental  measured  velocity  at  a point  and  the  computed  value. 
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Configuration 


Figure  36.  Faired  measured  axial  velocity  defect;  z = 3 in. 
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Figure  37-  Computed  sectional  profile  drag  coefficients  from 
Figures  35  and  36 
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Figure  30.  Configuration  1;  comparison  of  measured  and  computed 
horizontal  and  vertical  velocities  (x  * 50  in., 
z = 67.95  in.) 
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Figure  41.  Configuration  lj  measured  axial  velocities  (x  ■ 50  in.) 
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Figure  51.  Configuration  2;  flap  vortex.  Comparison  of  measured 
and  computed  horizontal  and  vertical  velocities 

( x » 50  in . , z = 6O.3  in . ) 


73 


. tpt 


z *61. 4 


-,ure  53-  Configuration  2;  measured  axial  velocities  (x  * 50  in.) 
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Figure  58.  Configuration  2;  flap  vortex.  Comparison  of  measured 
and  computed  horizontal  and  vertical  velocities 
(x  « me  in.,  z * 57.5  in. ) 
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Figure  61.  Configuration  2;  measured  axial  velocities  (x  = 100  in.) 


A weighting  function,  which  was  simply  the  radial  distance  from 
the  assumed  vortex  center  to  the  poirr  at  which  the  measurement 
was  made,  was  used  to  bias  the  most  distant  -a  (from  the  vortex 
center) . The  output  of  the  program  was  the  cortex  center  or 
centers  (when  both  a tip  and  flap  vortex  were  present).  When  the 
weighting  function  was  taken  to  be  unity,  the  vortex  centers  only 
changed  slightly.  The  horizontal  and  vertical  velocity  components 
were  then  computed  from  the  Betz  model  and  extension  for  flap 
vortices.  The  results  are  shown  in  Figures  38-61  and  in  Figures 
B.l  and  B.3.  The  calculation  was  performed  with  the  load  distri- 
butions obtained  from  pressure  measurements  (labeled  experimental) 
and  lifting  surface  theory  (labeled  theoretical). 


Computations  of  the  horizontal  and  vertical  velocity  compo- 
nents for  the  flap  vortices  of  configurations  3 and  4 were  not  in 
good  agreement  with  experimental  velocity  measurements  and  are  not 
shown  in  Figures  B.2  and  B.4.  The  explanation  of  this  discrepancy 
is  most  easily  seen  in  Figures  62  and  63.  Here  the  measured 
velocity  vectors  in  the  Trefftz  plane  located  at  6 chords  dewn- 
stream  have  been  plotted.  The  flap  vortex  velocity  tields  are 
markedly  nonsymmetric ; roll-up  apparently  is  not  complete.  Yates 
(Ref.  14)  has  shown  that  the  initial  rate  of  roll-up  is  strongly 
dependent  on  the  maximum  value  of  dr/dy  . Referring  back  to  the 
wing  lift  distributions  shown  in  Figures  24  and  25,  it  is  clear 
that  the  flap  vortices  will  develop  at  a lesser  rate  than  the  lip 
vortices.  Also,  it  is  expected  that  the  flap  vortex  in  configura- 
tion 2 will  develop  more  rapidly  than  those  of  configurations  3 
and  4.  In  Appendix  A,  a simple  model  which  estimates  rate  of  roll- 
up has  been  developed  for  simple  load  distributions.  This  mo. 'el 
also  shows  how  roll-up  rate  is  related  to  rate  of  change  of  the 
wing  load  distributions.  While  the  results  of  this  model,  in  its 
present  form,  cannot  be  directly  applied  to  flap  vortices,  it  is 
reasonable  to  assume  that  at  axial  stations  near  the  wing,  tip 
vortices  will,  in  general,  be  better  defined  (roll-up  will  be  more 
nearly  complete)  than  flap  vortices.  Figures  38—6 1 support  this 
assertion. 

4.  ROLL-UP  OF  THE  WING  DRAG  DISTRIBUTION 


The  vortex  axial  and  swirl  velocity  distributions  were 
calculated  for  the  clean  wing  configuration  omy,  U3lng  the  exten- 
ded Betz  model  developed  in  Section  11,2.  The  sectional  induced 
drag  coefficient  for  Cl  = 0.58  was  obtained  from  a lifting  surface 
calculation  and  is  shown  in  Figure  64.  The  section  drag  coeffi- 
cient is  the  sum  of  the  profile  drag  and  induced  drag  and  is  also 
shown  in  Figure  64.  The  results  of  the  calculation  for  the  axial 
and  swirl  velocity  distributions  are  shown  in  Figure  65.  As  can 
be  seen,  the  axial  velocity  is  nearly  the  free  stream  value.  This 
result  is  surprising  and  significant  in  that  it  3hows  that  the  drag 
on  a typical  airfoil  Is  nearly  that  required  to  keep  the  axial 
velocity  uniform  and  equal  to  the  free  stream  value  across  the 
vortex.  In  the  absence  of  drag,  the  axial  velocity  is  in  excess  of 
the  free  stream  value  as  discussed  by  Batchelor  (Ref.  17).  This 
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Measured  velocities  in  the  Trefftz  plane;  configuration  3 )z  * 50  in. 


Figure  63.  Measured  velocities  in  the  Trefftz  plane;  configurati 
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Figure  64.  The  sectional  Induced  drag  coefficient  for  configuration  1 (from  a lifting 
surface  calculation) 
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suggests  to  us  that  unless  large  drag-producing  devices,  such  as 
spoilers,  are  present  on  the  wing,  an  adequate  description  of  the 
downstream  vortex  may  be  obtained  by  setting  U(r,  * U«  in  Eqs. 
(13)  and  (14).  The  vortex  structure  is  then  given  by  the  Betz 
model  and  the  extension  given  by  Donaldson,  et  al.  There  is  often 
a strong  axial  velocity  excess  or  defect  on  the  vortex  centerline; 
however,  this  is  not  significant  in  that  the  actual  mass  flux  or 
axial  momentum  flux  departure  over  the  uniform  axial  velocity 
vortex  is  small  when  compared  to  these  quantities  over  the  whole 
vortex  (defined  as  that  region  in  which  all  the  axial  vortlcity 
may  be  found).  As  will  be  shown,  this  large  excess  or  deficit  on 
the  vortex  centerline  is  rapidly  diminished  by  viscous  processes. 

5.  ROLL-UP  OP  WING  TURBULENT  KINETIC  ENERGY 

The  initial  turbulent  kinetic  energy  distribution  in  the 
rolled-up  vortex  for  the  clean  wing  configuration  was  estimated 
with  the  procedure  given  in  Section  II. 3.  The  turbulent  kinetic 
energy  distribution  at  the  trailing  edge  of  the  wing  was  calcula- 
ted from  a turbulent  boundary  layer  program  developed  at  A.R.A.P. 
(Ref.  19).  The  distribution  of  q2  computed  at  the  trailing  edge 
is  shown  in  Figure  66a;  the  maximum  value  of  q/Uo>  was  nearly  0.09  • 
It  was  assumed  that  the  distribution  of  kinetic  energy  was  even  in 
z (where  z was  measured  vertically  from  the  wing).  The  distri- 
bution of  q’2  in  the  vortex  is  shown  in  Figure  66b.  Note  that 
the  maximum  (at  r » 0)  is  the  same  maximum  value  in  the  boundary 
layer.  The  sudden  drop-off,  as  mentioned  earlier,  is  brought  about 
by  the  mixing  of  turbulent  boundary  layer  fluid  with  increasing 
1 amounts  of  fluid  which  do  not  contain  turbulent  kinetic  energy  as 

loll-up  proceeds  from  the  tip. 

j 6.  DECAY  OF  AN  ISOLATED  TURBULENT  VORTEX 

A three-dimensional  turbulent  vortex  program  was  developed  at 
A.R.A.P.  under  contract  to  ARL.  This  program  is  described  in 
detail  in  Ref.  15.  Given  the  axial,  swirl,  and  second-order  velo- 
city correlation  radial  distributions,  the  program  will  calculate 
the  entire  mean  and  turbulent  vortex  structure  at  subsequent  posi- 
tions downstream.  Output  at  designated  downstream  stations  are 
the  mean  velocity  distributions  and  second-order  velocity  correla- 
tion distributions.  In  addition,  the  torque  exerted  by  the  vortex 
on  a flat  rectangular  airfoil  of  prescribed  span  may  be  calculated 
at  designated  downstream  distances.  We  have  run  the  program  for 
the  clean  wing  configuration  using  the  axial  and  swirl  velocities 
shown  in  Figure  65  as  calculated  from  the  experimentally  measured 
wing  load  distribution.  The  initial  turbulent  quantities  were 
specified  with  auto-correlatlons  taken  to  be  q,J(r)/3  (q,2(r) 

is  given  in  Figure  66)  and  cross-correlations  equal  to  zero. 

The  results  of  the  calculation  are  summarized  in  Figures  67- 
72.  Note  that  the  initial  swirl  and  axial  velocity  distributions 
were  smoothed  near  r = 0 to  avoid  numerical  problems  with  the 
number  of  grid  points  needed  to  adequately  represent  this  region. 
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Figure  66a.  Estimated  distribution  Oi  turbulent  kinetic  energy 
at  the  wing  trailing  edge 
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F.’gure  66b.  The  computed  distribution  of  turbulent  kinetic  energy 
in  the  Betz  vortex 
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In  general,  the  calculation  is  not  sensitive  to  this  type  of 
fairing  since  the  actual  axial  flux  of  linear  and  angular  momentum 
is  negligible  through  the  region. 

Two  results  of  this  calculation  are  significant.  The  first 
concerns  the  level  of  turbulence  which  is  supported  by  the  vortex. 
In  Figure  71  we  see  that  the  axial  flux  of  q'2  , after  the  adjust- 
ment of  the  inviscid  profiles  (two  semispans),  drops  far  below  the 
initial  value.  This  result  is  not  surprising  since,  if  the  produc- 
tion terms  for  the  second-order  correlations  are  examined  for  a 
turbulent  boundary  layer,  production  is  proportional  to  the  shear 
and  of  the  order  of  the  free  stream  velocity  divided  by  the 
boundary  layer  thickness.  In  a vortex  in  which  radial  gradients 
of  axial  velocity  are  small,  turbulent  production  is  proportional 
to  r(V/r)r  , a term  which  vanishes  in  the  viscous  core  of  the 
vortex  (where  V goes  like  r ) and  decays  as  r-2  in  the  irro- 
tational  region.  At  the  maximum  in  swirl  velocity,  production  is 
proportional  to  V / r , a quantity  which  is  small  when 

compared  to  on  the  wing.  The  fact  that  no  turbulent 

production  is  possible  at  the  vortex  centerline  suggests  an  expla- 
nation for  the  fact  that  the  maximum  value  of  q*  does  not  occur 
at  r = 0 as  shown  in  Figure  70. 

The  second  result  concerns  the  slow  decay  of  the  vortex  as 
measured  by  the  torque  calculated  on  a flat  plate  airfoil  (Fig. 

72).  We  do  not  want  this  result  to  imply  that  the  rate  of  decay 
of  an  isolated  vortex  represents  the  decay  rates  to  be  expected  in 
an  aircraft  wake.  The  assumption  of  axial  symmetry  and  assumptions 
regarding  the  nature  of  the  boundary  conditions  on  the  correlations 
for  large  r are  idealizations  which  become  invalid  as  the  vortex 
pair  ages.  We  note  that  the  calculation  was  carried  downstream  to 
a nondirnensional  distance  of  4o  semispans.  Investigations  to 
determine  the  amplification  rats  of  sinusoidal  instability  have 
shown  that  initial  disturbances  grow  to  e times  their  Initial 
amplitude  in  time  0'(27ib2/T)  or,  in  terms  of  downstream  distance, 
0'(8AR/Cjj)  wing  spans.  For  this  configuration,  the  downstream 
e-folding  distance  is  of  the  order  of  40  semispans,  and  at  40  semi- 
spans the  wake  is,  indeed,  it  ill  very  young. 

The  axial  velocity  distribution  shown  in  Figure  68  can  be 
compared  with  that  measured  for  the  clean  wing  (Fig.  45)  at  10 
chords  downstream.  The  comparison  is  made  iri  Figure  73  and,  as  can 
be  seen,  is  in  general  agreement. 


y,  inches 


Figure  73*  Comparison  of  the  computed  and  measured  axial 
velocities  at  x * 100  in.  (configuration  1) 


SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 


1.  CONCLUSIONS 

Several  conclusions  from  this  study  can  now  be  made  regarding 
the  structure  of  aircraft  wakes. 

1)  The  Betz  model  and  extension  to  include  "interior" 
vortices  give  a good  description  of  the  circulation 
distribution  in  wake  vortices  outside  the  viscous  core 
and  after  roll-up  is  complete.  At  distances  closer  to 
the  aircraft,  a more  detailed  model  is  needed.  One 
approach  would  be  to  treat  the  roll-up  region  as  a 
growing  vortex  and  stretching  sheet,  as  described  in 
Appendix  A.  At  distances  far  downstream,  aging  alters 
the  vortex  structure,  and  descriptions  of  the  decaying 
vortex  must  be  obtained  by  deliberate  calculation. 

2)  Distributed  wing  drag  of  the  magnitude  typical  of  clean 
airfoil  sections  has  only  a small  effect  on  the  inviscid 
tip  vortex  structure.  This  results  primarily  from  the 
fact  that  typical  values  of  drag  reduce  the  axial 
velocity  in  the  vortex  so  that  it  is  nearly  uniform 
and  equal  to  the  free  stream  value.  This  result,  while 
obtained  by  direct  calculation  for  a tip  vortex,  is 
expected  to  be  approximately  correct  for  "interior" 
roll-ups  as  well.  When  large  drag  devices  are  deployed 
on  a wing,  significant  structural  changes  in  the  wake 
vortex  can  result.  Drag  devices  which  leave  the  lift 
distribution  unchanged  result  in  vortex  structural 
changes  which  are  a consequerce  of  redistributing  the 
axial  flux  of  angular  momentum  in  the  vortex.  The 
total  flux  of  angular  momentum  remains  unchanged  in  the 
region  containing  axial  vorticity.  The  inviscia  model 
suggests  that  vortex  deintensification  broueht  about 
solely  by  increased  drag  will  be  expensive  in  terms  of 
drag  penalty.  It  is  likely  that  drag  devices  will 
significantly  raise  the  turbulent  level  in  the  wake 
vortex  which  will  lead  to  increased  rates  at  which  the  . 
angular  momentum  is  redistributed  outward.  Therefore, 
estimates  of  the  power  required  to  reduce  the  intensity 
of  a wake  vortex,  as  calculated  from  an  inviscid  model, 
are  likely  to  be  high. 

3)  Based  on  the  models  used  in  this  study,  the  flux  of 
turbulent  kinetic  energy  in  the  wing  boundary  layer 
which  rolls-up,  into  a vortex  is  greater  than  the  flux 
of  turbulent  kinetic  energy  that  is  to  be  found  in  the 
downstream  vortex.  Since  asymmetries  exist  in  the  flow 
field  of  a vortex  pair  or  pairs,  it  is  likely  that  aging 
will  occur  more  rapidly  than  is  predicted  here  by 
computing  the  decay  of  an  isolated  turbulent  vortex. 
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1} ) The  technique  developed  to  calculate  the  location  of 

the  centroids  of  the  discrete  concentrations  of  vorticity 
in  the  wake  have  been  shown  to  be  quite  accurate  by 
comparison  with  flow  visualization  studies  made  here  and 
by  others.  It  has  been  shown  that  the  structure  of  two- 
vortex-pair  wakes  can  be  obtained  from  a simple  classi- 
fication chart.  Descriptions  of  the  relative  positions 
of  the  centroids  of  three  of  more  pair  wakes  must  be 
obtained,  in  general,  by  direct  calculation. 

2.  RECOMMENDATIONS 

Our  recommendations  cf  additional  work  which  should  be  under- 
taken in  the  area  of  vortex  wake  dynamics  are  divided  into  two 
categories : 

1)  Application  of  the  technology  de\eloped  in  this  and 
previous  studies; 

2)  Additional  research  to  understand  the  complete  vortex 
wake  behavior. 

1)  Application  of  New  Technology 

It  has  been  shown  that  in  the  spirit  of  Betz  mean  discrete 
vortex  structure  may  be  calculated  from  the  wing  span  load  and  drag 
distributions,  and  a simple  model  may  be  used  to  calculate  the 
relative  positions  of  these  vortices.  This  In  Itself  provides  a 
description  of  the  rolled-up  aircraft  wake  with  a precision  that 
was  previously  unavailable.  In  addition,  the  model  developed  in 
Appendix  A provides  a first  step  In  a description  of  the  wake  during 
roll-up,  from  which  the  velocity  field  in  the  vicinity  of  the  wing 
may  now  be  calculated  in  detail.  These  models  have  Immediate  appli- 
cation to  the  design  and  evaluation  of  the  performance  and  stability 
of  lifting  surfaces  in  vortex  wake  flow  fields.  Some  obvious 
applications  of  these  wake  models  are 

1)  The  design  of  tails,  canards,  side  force  generators,  and 
w^ngs  In  vortex  flow  fields  produced  by  upstream  lifting 
surfaces.  The  longitudinal  instability  of  the  E.B.F.  V/ 

STOL  model  is  but  one  example  of  the  need  for  the 
application  of  this  new  technology. 

2)  The  use  of  velocity  fields  as  input  to  flight  simulators. 
Pilots  could  then  train  for  wake  encounters;  military 
pilots  could  simulate  aircraft  upset  due  to  wake  turbu- 
lence under  close-interval  take-off,  formation  flying, 
and  refueling  conditions. 

3)  The  design  of  a low  hazard  wake.  By  specifying  the 
vortex  wake  structure,  the  wing  lift  and  drag  distri- 
bution which  generates  the  desired  wake  may  be  calculated. 
Estimates  of  the  induced  and  profile  drag  penalties  to 
obtain  a specified  level  of  wake  intensity  may  be 
obtained  by  a direct  calculation. 
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2)  Additional  Research 

The  wake  models  presented  in  this  study  provide  an  important 
step  in  developing  techniques  to  describe  the  complete  vortex  wake 
history.  However,  several  important  aspects  of  the  problem  still 
remain  unsolved.  We  have  proposed  that  a complete  description  of 
the  aircraft  wake  must  answer  three  basic  questions  — roll-up, 
aging,  and  Instability.  The  Betz  model  and  extensions,  including 
the  roll-up  model,  can  only  be  used  with  confidence  before  aging 
and  Instability  become  important.  The  areas  of  wake  aging  and 
vortex  breakdown  need  additional  study. 

Our  understanding  of  the  phenomenon  of  vortex  pair  or  pairs 
aging  is  Just  in  its  infancy.  Although  fundamental  contributions 
have  been  made  to  our  understanding  of  the  aging  of  an  axisymmet- 
ric  isolated  turbulent  line  vortex,  few  of  the  conclusions  which 
may  be  drawn  can  be  directly  applied  to  the  pair.  Advances  in 
numerical  computation  procedures  ft  A.R.A.P.  and  increased  confi- 
dence in  the  invariant  modeling  technique  have  recently  made  the 
calculation  of  turbulent  decay  of  a complete  vortex  wake  possible. 
We  strongly  recommend  that  this  unique  capability  be  exploited  to 
make  more  realistic  calculations  of  wake  decay  than  have  previous- 
ly been  possible. 

Vortex  breakdown  occurring  in  proximity  to  lifting  surfaces 
can  have  adverse  elfects  on  the  performance  of  those  surfaces. 

New  advances  in  our  understanding  of  the  phenomenon  can  now  predict 
conditions  conducive  to  breakdown.  As  advanced  tecnnology  aircraft 
use  favorable  lifting  surface  vortex  interaction  effects  to 
improve  performance,  the  ability  to  predict  and  control  breakuown 
will  have  iirfcortant  consequences.  Favorable  vortex-surface  inter- 
action, when  exploited,  will  undoubtedly  improve  many  aspects  of 
aircraft  performance , 
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APPENDIX  A 


ESTIMATED  TIME  TO  ROLL  UP  A TWO-DIMENSIONAL  SHEET 


The  assumptions  made  by  Betz  can  be  used  to  develop  a model 
to  estimate  the  time  required  to  roll  up  the  vorticity  behind  a 
simply  loaded  wing.  In  the  spirit  of  Betz,  we  consider  the  roll- 
up of  a two-dimensional  vortex  sheet  of  strength  y(y»t). 


Referring  to  Figure  A-l,  we  assume  that  vorticity  which  has 
been  convected  past  station  A is  rolled  up  in  the  vortex  end  is 
distributed  according  to  Eq.  (^).  Further,  we  assume  that  the 
portion  of  the  sheet  which  is  not  yet  in  the  vortex  remains 
horizontal  as  shown.  The  error  which  is  a consequence  of  this 
assumption  is  discussed  below. 


The  time  rate  of  change  of  circulation  in  the  vortex  is 
given  by 


df 1 
dt 


Y(yv»t) 


( A . 1 ) 


where  v is  the  horizontal  velocity  at  station  A induced  by  the 
vorticity  and  y is  the  vortex  sheet  strength  at  station  A.  The 
term  dyy/dt  accounts  for  the  inward  motion  of  the  vortex  as 
roll-up  proceeds. 

In  the  evaluation  of  the  sheet  strength,  it  is  necessary  to 
account  for  sheet  stretching  which  is  the  result  of  nonuniform 
convection  along  the  sheet.  To  obtain  an  expression  governing  the 
stretching,  the  two-dimensional  inviscid  vorticity  equatjon 
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is  integrated  through  the  sheet  from  z = -»  to  z * ».  After 
integrating  by  parts  and  using  the  continuity  equation,  we  obtain 
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(A. 3) 


where  y = / u>  dz  . 

-00 

The  horizontal  velocity  v may  be  thought  of  as  the  sum  of 
two  terms  vc  and  vs  ; vc  is  the  velocity  due  to  the  rolled-up 
portion  of  the  sheet  and  is  given  by 

v0  ■ §r  { -[<y  - ?v)2  + z2]  + [(y  + yv)2  + z2]  (A-1|) 
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When  a vortex  sheet  is  planar,  vs  is  identically  zero  which 
motivated  the  assumption  regarding  the  geometry  of  the  sheet.  The 
error  introduced  is  not  too  severe,  however,  since  had  the  sheet 
been  allowed  to  deform  having  a characteristic  radius  R , the 
horizontal  induced  velocities  would  be  of  order  Fs/R  where  rs 
is  the  circulation  of  the  sheet.  Initially,  R is  infinite.  When 
roll-up  has  proceeded  for  some  time,  R is  finite;  however,  Ts 
is  small  since  most  of  the  circulation  has  already  been  rolled  up 
into  the  vortex.  We  therefore  neglect  the  right-hand  side  of  (A. 5). 
This  approximation  implies  that  the  roll-up  phenomenon  is  strongly 
dominated  by  the  developing  vortex.  The  approximations  Introduced 
thus  far  will  underestimate  the  time  to  roll-up,  since  sheet 
stretching  is  also  underestimate!. 


Unfortunately,  even  with  the  above-noted  simplifications,  the 
solution  of  (A. 5)  with  vc  given  by  (A.*0  is  quite  difficult. 
Therefore,  the  functional  form  of  the  convecting  velocity  is  chosen 
so  that  it  is  possible  to  find  an  analytic  solution  to  (A. 5)  and 
yet  retain  the  physics  of  the  stretching  phenomenon.  One  such 
velocity  is 

vc  = g(t)y 

where 

g(t)  ' 


(A. 6) 
(A. 7) 


The  general  solution  of  (A. 5)  can  now  be  found  by  the  method  of 
characteristics  and  is 


Y(y,t) 


rof<6> 
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(A. 8) 


where 


exp(- 


g(t)dt) 


(A. 9) 


f is  an  arbitrary  function  of  6 and  is  determined  from  the 
initial  sheet  strength  distribution. 


We  will  present  results  for  the  linear,  parabolic,  and  ellip- 
tical wing  loading.  The  initial  sheet  strengths  are  calculated 
from  y * -df/dy  and  are 
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[ Elliptic 


The  vortex  sheet  strength  at  station  A as  a function  of  time  Is 

f E(t) 
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The  location  of  the  developing  vortex  yv(l ) is  not  known  and 
Is  determined  from  the  conservation  of  impulse.  The  calculation 
appears  tractable  for  only  the  linear  load  case  where 


yv(t) 
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(A. 13) 


However,  a good  approximation  to  yv  is  to  use  y as  defined  in 
Eq.  (2)  or  rewritten  as  an  explicit  function  of  f as 


y(r)  = y(D 

For  the  linear  load,  y/s  * 1 - r/2f 


- * 


(A. 14) 
(A. 15) 


When  (A. 15)  is  used  to  approximate  yv  , an  error  of  6.25 % in  the 
impulse  for  the  linear  load  occurs  during  the  roll-up.  The  error 
is  less  than  this  value  for  lift  distributions  more  highly  loaded 
at  the  tip  and  is,  therefore,  not  serious. 


For  wings  having  a linear  load  distribution,  an  ordinary  non- 
linear integrodif ferential  equation  describing  the  time  rate  of 
growth  of  the  vortex  is 


dP 

dt 


l-l  o 
dt 
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(A. 16) 
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The  solution  is 


(A. 17) 


The  equations  governing  the  roll-up  for  the  parabolic  and  elliptic- 
ally  loaded  wing  are  considerably  more  complicated  and  are  not 
given  here.  The  solutions  for  these  cases  must  be  determined 
numerically  and  are  shown  along  with  the  linear  case  in  Figure  A-2. 

Spreiter  and  Sacks  (Ref.  7)  have  also  made  estimates  of  the 
downstream  distance  at  which  the  vortex  can  be  considered  to  be 
essentially  rolled  up  by  applying  the  results  obtained  by  Kaden 
(Ref.  6)  for  a semi-infinite  wing.  For  elliptic  wing  loading,  they 
obtained  the  nondimens ional  downstream  distance  of  xCjy^lyAR  * 0.18 
which  corresponds  to  having  approximately  55%  of  t.he  wing  root 
circulation  in  the  vortex  (as  calculated  here).  The  discrepancy  is 
believed  to  arise  as  a consequence  of  using  Kaden' s solution  where 
it  is  not  strictly  valid. 

For  small  times,  the  solutions  can  be  shown  to  behave  as 
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The  t J behavior  for  the  elliptic  load  agrees  with  the  similar- 
ity solution  obtained  by  Kaden. 

In  Figure  A-3  our  results  are  compared  with  those  recently 
obtained  by  Moore  (Ref.  5).  Part  of  the  discrepancy  has  already 
been  explained  by  the  approximations  vhich  have  been  made.  However, 
an  additional  point  is  that  Moore  has  chosen  the  station  past  which 
the  sheet  is  to  be  considered  rolled  up  at  a location  90°  in  the 
counterclockwise  sense  from  station  A in  Figure  A-l.  His  results 
are,  therefore,  biased  to  be  lower  than  those  obtained  here.  The 
amount  Is  difficult  to  calculate,  but  the  difference  is  expected  to 
be  most  significant  for  small  times.  The  result,  neglecting  sheet 
stretching,  is  obtained  by  taking  vc  constant  in  Eq,  (A. 5)  equal 
to  the  horizontal  velocity  at  station  A in  Figure  A-l,  In  this 
case,  complete  roll-up  occurs  in  finite  time. 
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Figure  A-2 . Vortex  circulation  and  radios  as  a function  of  downstream  distance 
(time  has  been  replaced  by  x/U*) 


APPENDIX  B 


ADDITIONAL  EXPERIMENTAL  WAKE  VELOCITY  DISTRIBUTIONS 


Appendix  B contains  measured  wake  vortex  velocity  distributions 
for  configurations  3 and  4.  Swirl  velocity  distributions  are 
computed  and  compared  for  tip  vortices  only.  Flap  vortices  have 
not  rolled  up  sufficiently  at  the  five-chord  downstream  measuring 
station  to  discern  a distinct  axisymmetric  structure. 
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Figure  3-la.  Configuration  3;  tip  vortex.  Comparison  of  measured 
and  computed  horizontal  and  vertical  velocities 
(x  = 50  in. , z - 67 in. ) 
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Figure 


B-2b , Configuration  3;  flap  vortex.  Measured  horizontal 
and  vertical  velocities  (x  = 50  in.,  z » 67.4  in.) 
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Figure  B-2d.  Configuration  3;  measured  axial  velocities  (x  * 50  in.) 
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Figure 


-3a.  Configuration  tip  vortex.  Comparison  of  measured 
and  computed  horizontal  and  vertical  velocities 
(x  = 50  in.,  z * 5^.5  in.) 
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Figure  S3-3d.  Configuration  measured  axial  velocities  (x=  50  in.) 

123 


r Figure  B-4b.  Configuration  4;  flap  vortex.  Measured  horizontal 

1 and  vertical  velocities  (x  “ 50  in.,  z = 49.9  in.) 
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